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ABSTRACT 

We present initial results from time series imaging at infrared wavelengths of 
0.9 sq. degrees in the Orion Nebula Cluster (ONC). During Fall 2009 we obtained 
81 epochs of Spitzer 3.6 and 4.5 /um data over 40 consecutive days. We extracted 
light curves with ~3% photometric accuracy for ~2000 ONC members ranging 
from several solar masses down to well below the hydrogen burning mass limit. 
For many of the stars, we also have time-series photometry obtained at optical (Jc) 
and/or near-infrared (JK^) wavelengths. Our data set can be mined to determine 
stellar rotation periods, identify new pre-main-sequence (PMS) eclipsing binaries, 
search for new substellar Orion members, and help better determine the frequency 
of circumstellar disks as a function of stellar mass in the ONC. Our primary 
focus is the unique ability of 3.6 & 4.5 //m variability information to improve our 
understanding of inner disk processes and structure in the Class I and II young 
stellar objects (YSOs). In this paper, we provide a brief overview of the YSOVAR 
Orion data obtained in Fall 2009, and we highlight our light curves for AA-Tau 
analogs - YSOs with narrow dips in flux, most probably due to disk density 
structures passing through our line of sight. Detailed follow-up observations are 
needed in order to better quantify the nature of the obscuring bodies and what 
this implies for the structure of the inner disks of YSOs. 

Subject headings: open clusters and associations: individual (Orion) — circumstellar 
matter — stars: pre-main sequence — stars: protostars — stars: variables: general 



1. Introduction 



The Orion star forming region has bee n the s ubjec t of more low mass young star var i- 
abil ity studies tha n any other region. While [Hard (119691 ). iHerbig fc Kameswara Rad ( 1l972l ). 
and IWalken (119781 ) all conducted pioneering studies of the variability of young stars in Orion, 
the modern era usirig two - dimensional imaging cameras was initiated by Herbst and his team 
( lAttridge fc Herbstlll992l : IChoi fc Herbstlll996l ). They obtained multi-year, optical time se- 
ries photometry of several regions of the ONC, eventually obtaining light curves for hundreds 
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of PMS stars, often spanning several years. Those data were primarily used to derive rotation 
periods - made possible by the fact that the photospheres of young stars are often heavily 
spotted (cold and/or hot spots), resulting in flux variations modulated at the star's rotation 
period. Many other groups subsequently obtained time series photometry of other portions 



of th e ONC, using o ptical ([Stassun et al. 



999; RebuU 2001; Herbst et al. 2002; Irwin et al. 



20071 ) and near-IR ( Carpenter et al. booif T imaging data. In many cases, the light curve 



shapes are well fit by hot or cold spots, normally at moderately high latitudes since the light 
curves are seldom "fiat-bottomed." However, for a significant fraction of the light curves, 
particularly those in the near-IR, spots do not seem to provide a good explanation of the 
observed variability (ICarpenter et al.ll200l[ ). 



Time series photometry of YSO s can also address issue s other than rotatio nal velocities 



2007: Herbst et al. 



12010 



such as fiares and fiare f requencv (IParsamian et al.lll993l ). disk instabilities ( iFedele et al. 



Plavchan et al. 



2008ar). and echpsing binaries as tests of PMS 



isochrones ( IStassun et al.ll2004j . l2006l : ICargile et al.ll2008l ). Mid-IR time series photometry of 
YSOs could, in principle, detect variations in the temperature or spatial distribution of warm 
dust in the inner circumstellar disks of these stars, thereby informing models of disk evolu- 
tion and perhaps of planet formation. We recently completed the first sensitive , wide-area. 



mid-I R time-series photometric moni toring program for the ONC using IRAC (IFazio et al. 



2004) on the Spitzer Space Telescope (jWerner et al. 2004). We report here some early results 
from that program, with more detailed analysis of specific classes of variability to follow. 



2. Observations 

2.1. Spitzer Observations 

Our GOG Exploration Science program YSOVAR (Young Stellar Object VARiability) 
provides the first large-scale survey of mid-IR photometric variability of YSOs. In Fall 2009, 
we used ~250 hours of warm Spitzer observing time to monitor ~0.9 deg^ of the ONC at 
3.6 and 4.5 /im. 

The observed area was broken into five segments with a central region of ~20' x 25' and 
four fianking fields. The central part was observed in full array mode with 1.2 seconds expo- 
sure time and 20 dithering positions to avoid saturation by the bright nebulosity around the 
Trapezium stars. The remaining four segments of the map were observed in High Dynamic 
Range mode with exposure times of 10.4 and 0.4 seconds, and 4 dithering positions. A 
summary of the observations appears in Table [U 

These observations were taken for 40 days in Fall 2009, with ~2 epochs each day but 
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with an interval between observations that varied both by design (to reduce period-ahasing 
problems) and due to constraints imposed by other Spitzer programs th at were executed 



durin g the same campaigns. We used the IDL package Cluster Grinder (iGutermuth et al. 



20091 ) which, starting from the BCD images released by the Spitzer Science Center, builds 



the combined mosaic for each epoch and performs aperture photometry on the mosaics. In 
the present work the BCD images used correspond to the S18.12 software build (the zero 
point magnitudes are 19.30 and 18.67 for 1 DN/s total flux at 3.6 and 4.5 /im, respectively). 
Future versions of the pipeline will include improvements, most importantly using better 
linearity corrections. 



2.2. Ground-based Data 

In order to complement our Spitzer data, we also obtained Ic, J, and Kg contemporane- 
ous photometry, usually for smaller areas within the Spitzer mosaic. The main datasets come 
from the UKIRT Wide Field CAMera (WFCAM), the CFHT Wide field InfraRed Camera 
(WIRCam), the Steward Observatory Super-LOTIS (Livermore Optical Transient Imaging 
System) robotic telescope, and the NMSU/APO Im telescope. A summary of the principal 
characteristics of each dataset appears in Table [H A deeper analysis on data reduction of 
this data set plus additional data collected for this project will be described in Gutermuth 
et al. (2011, in prep.). 

We performed differential aperture photometry on the ground-based data. In each 
dataset, for each object, an artificial reference level was created by selecting 30 nearby 
isolated stars; we iteratively eliminated those with larger photometric scatter or evidence 
of variability. Finally, the time series are incorporated into a database that will become 
publicly available in early 2011 (http:/ /ysovar.ipac.c altech.edu/) . This database includes: 
a color image of the surveyed field of view, fits files for the 3.6 and 4.5 ^m mosaics from a 
single epoch, a table for all ONC stars providing the name, aliases, RA, Dec, J,H,K, [3.6], 
[4.5], [5.8], [8.0], and [24] information as in Megeath et al. (2011, in prep.), a file with plots 
of multi-wavelength light curves, and the actual time series for all ONC members either as 
ascii files or by querying the database. Alternatively, the time series can be found in Tabl^ 



3. Variability statistics and types 



Light curves were extracted in both IRAC bands for 1249 Class I and II previously known 
Orion YSOs with IR excess (Megeath et al. 2011, in prep.) and inferred masses between 
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0.03 and ~2 M0, plus 820 other likely Orion members. This latter set of objects is formed by 
sources that fulfill at least one of the following c o nditions: a) has been labeled as a member b y 



proper mot i on stu dies such as iParenagd (Il954j ) ; iJones fc Walker! (I19881): iMcNamara 



Tian et al.l (Il996l ). b) has been labeled as a membe r according to 



c) has been labeled as a member by X-r av studies (|Getmari et al. 



typical of weak-lined T Tauris (WTTs ) (ITobin et al. 



Hillenbrand et al. 



2009 



1976 ): 



f ll998[ ) 



2005h d) h as H„ profiles 
Fiiresz et al.l 120081 ) . e) has been 



labeled as variable by ICarpenter et al.l ( 1200 ll ). or f) finally we have included all the sources 
that have been claimed to be periodic variables. Most of these stars are WTTs according 
to the Megeath et al. photometry, though some may be sources with excesses that escaped 
previous detection. In addition to these 2069 Orion stars, another 147 Class I and Class II 
sources plus 209 Orion members without IR excess have light curves in just one of the two 
IRAC bands. 

We estimated the noise in the IRAC light curves by measuring the typical scatter for all 
the objects in the field. In that way, we estimate a relative per channel photometric accuracy 
better than 3% down to [3.6]=14 and better than 10% for objects as faint as [3.6]=16. For 
the J band, we derive accuracies of ~2% down to J=16.5 and of ~10% to J=18.5 for the 
ensemble. Note that the accuracy varies significantly over the surveyed field due to nebular 
emission. 

Approximately 70% of the YSOs with IR excess are var iable based primarily on the 
Stetson index calculated using both IRAC channels together (IWelch &: StetsonI (119931 ): see 
Table m for statistics). The variable sources were visually inspected and, in addition to the 
Stetson variables, we considered as variables as well a few YSOs for which a period could be 
found even when the Stetson index alone would not reflect the variability. About 3% of the 
variable YSOs have been selected in this way. Usually these sources have very low amplitude 
however, most of them have previously detected periods that match the ones derived here and 
thus we believe the variability is real. The variable stars exhibit diverse behavior with time, 
including slow drifts in brightness over one month or longer that may be caused by a slow 
change in accretion rate or by self-shadowing in a warped disk (Fig. [1^); rapid chaotic flux 
variations probably caused by changes in the accretion geometry (Fig. [T]d); flares (Fig. [It); 
stars with lightcurves that look "periodic" for part but not all of the 40 day period or for 
which the periodic variations are mixed with other longer term changes (Fig. [T]i); periodic 
or pseudo-periodic variations arising from photospheric spots and/or disk warps (Fig. [1^); 
changes in color from variable extinction (?) (Fig. [T]F); and stars that show narrow or broad 
dips in the observed flux superimposed on an otherwise nearly constant flux level (see next 
section) . 



Fig. [2] illustrates the value of having contemporaneous optical/NIR/Mid-IR monitoring 
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(12%, 85%, and 37% of the Orion stars have K^, J and Ic band hght curves respectively). 
When compared to the near-IR and optical data, the mid-IR light curves usually have similar 
shapes. We have derived the amplitude of the light curves at each band as the difference be- 
tween the maximum and minimum values of the magnitude at that bandpass avoiding single 
deviant datapoints. There are generally larger amplitudes at shorter wavelengths (~60% of 
the variables with shorter wavelength data- Fig. [2^), though in some cases, the amplitudes 
of variation are similar at all bands (~30% of the variables with shorter wavelength data- 
Fig. Wp)- In a small number of cases, the IRAC variations are smaller or non-apparent at 
shorter wavelengths, indicating that the source of variability is probably in the disk (~10% 
of the variables with shorter wavelength data- Fig. [S^d). In only two cases do we find a 
phase difference between the optical and near-IR data relative to IRAC data (Fig. [2^). Some 
of these effects can be explained by the presence of one or more hot spots on the surface of 
the star and/or the existence of warps in its disk. 



Several of the objects deserve particular mention. Fig. [T]i shows the we 



1-known inter- 



mediate mass protostar OMC-2 IRS 1, f or wh i ch outflows have been detected (|Yu et al.lll997l : 
Reipurth etaU Il999l : IXakahashi et all boosi : lAnathpindika fc WhitworthI boosh: a similar 
l ight c urve was previously reported for Orion BN/KL fJ053514. 12-052222.9 - iHillenbrand et al 



( 2001)), though the physical mechanism for such variability in these B-type (jjohnson et al. 



1990l : IHillenbrand et al.ll200ll ) YSOs is not obvio us. In Fig . [2b, we show the protostar OMC-2 
IRS2, which has also been associated with a jet (lYu et al.lll997l : lAnathpindika &: Whitworth 
20081 ) but for which only one sca ttering lobe is seen , suggesting that its disk is significantly 
less inclined to the line-of- s ight (iRayner et al.lll989[ ). Fig. |2]F shows the lightcurve of a star 
identified by IStassun et al.l (119991 ) as the fastest rotator in a s tudy of the period distribution 
in the ONC (P=0.27 days. Star 1161 of Istassun et"aD Jl999 ) (RA 05 34 45.93, Dec -04 49 
22.0, from his Table 1). Given its spectral type (K5; iRebullI (120011 )) . at 1-2 Myr, this period 
is shorter than break-up. However, the VIcJHK photometry for the star suggest an age 
ol der than 30 M yr if it is at Orion distance. Because the light curve shape and amplitude 



m 



Stassun et al.l appears identical to our 2009 light curve, we suspect this object is a tidally 



distorted binary, but it could also be a single near ZAMS rapid rotator (similar to HII1883 
in the Pleiades). 

In general, our variable YSOs show colorless changes at IRAC bands (~88% of our 
sample show IRAC [3.6] — [4.5] RMS amplitudes smaller than 0.05 mag). The objects that 
do show IRAC color variations usually become redder as they fade. The typical peak-to-peak 
amplitudes in each IRAC band are ~0.2 mag. The most extreme stars in our sample have 
changes in amplitude as large as ~1.8 mag and the largest change in color is ~0.3 mag. 

We have used the NStED Periodogram ServiceBto derive periods for our sample. This 
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service uses the Lomb-Scargle (Scarglelll982[) . Box-fitting Least Squares (IKovacs et al.ll2002[ ) 



and Plavchan ( Plavchan et al. 2008bl ) algorithms for computing periodograms from time 



series data. Given the cadence of our data, the periodograms are formally reliable for periods 
between 1 and 40 days however we have found 11 sources with periods shorter than a day. 
Those sources all have periods from the literature that match our determinations and the 
light curves are so well phased that we believe the periods are real. We have not retained 
any period longer than 25 days. In this way we can determine a period for 23% of the 
variable YSOs with IR excess; 7% of these are Class Is and 93% are Class lis. Among the 
non-excess sources, ~44% are variables and the predominant variablity is periodic, probably 
due to photospheric spots; 33% of the WTTs have periods discernible from the Spitzer 
data. In the total ONC sample (1249 Spitzer excess sources plus 820 Orion members not 
known to have IR excess), we find 150 new periods (see Tabl^. About 30% of our total 
sample had previously published periods. Of these, we recover the published period for 57% 
of the WTTs and 47% for the YSOs with IR excess. This difference is likely due to the 
thermal emission from dust at IRAC wavelengths significantly degrading th e sensitivity to 
photospheric phenomena but also because the stars themselves may change ( lRebullll200ll ). 



Finally, we ha ve identified 9 eclip sing binarie s in our sample. Fo ur of ther n, 2MASS 

20061) ■ VI 174 Ori Jstassun et al.ll2004l ). Parl802 dCargile et al" 



05352 184-0546085 fIStassun et a 



20081), and JW 380 f llrwin et aL 



2007D are already known. A fifth one, Theta 1 Ori E, is a 



known spectroscopic binary (IHerbig fc GrifSnl 120061 ) which was flagged as potentially eclips- 
ing, but no confirmatory photometry had been obtained until now. We are working to obtain 
spectroscopic confirmation of the four new eclipsing binary candidates and will report on 
their status in a future paper. Their coordinates and tentative periods can be found in 
Table |5l 



4. "Dipper" Objects: Clouds in the disk? 



4.1. Identification and Characterization 



One of the most surprising and interesting classes of variability we find is narrow flux 
dips. This subset of 41 variables shows brief, sharp drops in their flux density with durations 
of one to a few d ays. T hese objects are members of the "Type III" variable class identified 
bv iHerbst et al.l (|l994f) based on an optical synoptic data set; their most famous prototype 
is AA Tau teouvier et al. 2003). A few examples of such "dipper" light curves are shown 



*http:/ /nsted. ipac.caltech.edu/periodograni/cgi-bin/Periodograin/nph-sinipleupload 
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in Figure 131 Among the identified cases, some exhibit only one dip, but typically >1 dip is 
seen over the 40 days. Periodic dips are also detected in about a third of these objects, with 
periods from ~2 to 14 days. The dip amplitude is typically several tenths of a magnitude 
in the IRAC bands and ranges from 0.16-1.5 mag at J. 

The criteria imposed for selecting dipper stars were: 1) a distinctive dip extends over 
several IRAC epochs unless a single epoch IRAC dip is corroborated by the independent 
but contemporaneous shorter wavelength data. Implicitly, this means the IRAC-only dips 
will be broader in time, on average, than dips where shorter wavelength data confirm the 
dip. 2) If based on just an apparent dip in one IRAC epoch, both IRAC channels must 
agree and the Ic or J-band data must show a larger amplitude dip by at least 50% to avoid 
misidentification of eclipsing binary stars as AA Tau analogs. 3) The "continuum" must 
be flat enough so that the IRAC dip "stands out." Although their identification by eye is 
generally straightforward, it is not as straightforward to determine the frequency of such 
dips in the overall data set. To do so, one must account for the likely superposition of such 
events with other types of variability that are occurring in the same stars. The optical and 
near-IR data were obtained at intervals of one to a few days as compared with the IRAC 
data which were obtained at half a day intervals, thus preventing confirmation of some IRAC 
dips. To detect these events requires good S/N data in all the relevant bandpasses, which 
eliminates many of the fainter stars and/or stars in regions where the nebular background 
is bright or very structured. 

Based on the above criteria, we identify 38 of the Class I or II YSOs as having one or 
more AA Tau-like flux dips. Given the coverage, cadence and S/N that we have in different 
filters, we estimate that we would be able to identify ~70% of broad dips (Fig. |3Kb) and 
about 40% of the narrow dips (those present in just one IRAC epoch; Fig. [3lie). Thus, at 
present, we estimate the overall frequency of dipp er stars at ^5% or hig her. We note that this 



fraction is much lower than the 28% estimated by lAlencar et al.l (120101 ) for periodic "AA Tau 
like" behavior in optical wavelength CoRoT data among a smaller set of young stars in NGC 
2264. There are, in fact, a variety of reasons that CoRoT detects a larger fraction of "dipper" 
sources in NGC 2264 than the Spitzer YSOVAR survey of the ONC: CoRoT's much better 
relative photometric accuracy (< 1 mmag) and higher cadence (every 512 sec during 23 days 
of uninterrupted observations), the larger amplitude of these dips at shorter wavelengths, the 
selection criteria itself, and the greater complexity of the mid-IR light curves due to the disk 
contribution. However, ther e may also be some rea l difference between the YSOs in the two 
clusters. Examination of the lCarpenter et al.l (120011 ) monitoring NIR data reveals that ~50% 



of our dippers displayed similar behavior 10 years ago, suggesting some stability in physical 
mechanism inducing the variations. We have identified 13 objects, or 34% of the "dipper 



sample, as being periodic with P<40 days (the observing window). ICarpenter et al.l (120011 ) 
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derived similar periods (less than 0.6 days of difference and usually less than O.ldays) for five 
of these objects. The dependence of depth vs. wavelength is roughly consis tent with that 



expec ted for the standard extinction law, as previously reported for AA Tau (IBouvier et al. 



20031 ). Coordinates, periods, and amplitudes for the dippers can be found in Table O 



Considering other existing data for our ONC stars, specifically the Spitzer photometry 
including all IRAC bands and MIPS 24 /im, 2MASS near-infrared, and literature optical 
photometry and spectral types, we find that there is no clear distinction of the AA Tau 
analogs in either Spitzer/IRAC or 2MASS magnitudes and colors relative to the rest of the 
ONC population, other than these objects are most often associated with Class II type YSOs. 
Only one of the 129 Class Is vs. 37 of the 1175 Class lis are identified as AA-Tau analogs. 



4.2. Interpretation 



The AA Tau variability phenomenon has been interpreted (jBertoutll2000l : iBouvier et al. 



20031 ) as arising from the rotation of a circumstellar disk with a high latitude "warp" that 
periodically obscures the star in addition to the extinction due to a flared disk that is 
typically imposed over the rest of the orbit. However, any process which creates over-dense 
asymmetric regions in the inner disk could produce the flux dips. 

Our interpretation is that these objects are being extincted by either "clouds" of rela- 
tively higher opacity in the disk atmosphere or geometric warps in the inner disk. In either 
case the dips are caused as the feature passes through our hne of sight to the star as the disk 
rotates. This hypothesis is consistent with both the event duration and the variation of dip 
amplitude with wavelength. Further, as some of the dips are periodic or semi-periodic, we 
note that the derived periods are consistent with those expected from an inner disk region in 
co-rotation with a star having typical rotation rates for a Class II star (2-14 days). Similar 
features located further out in the disk could be responsible for the non-repeating and/or 
broader dips. 

As an example, for a period of P=8 days, the "cloud" , a region of enhanced dust column 
density, shou ld be located at a K eplerian orbital radius 0.07 AU (or 7 R*, for a typical 1 



Myr M2 star; iBaraffe et al.l ( 119981 )) and for P>40 days (i.e., objects showing just one event), 
the occulting source would be located further than 0.19 AU (or 18 R,,). We can roughly 
estimate the size and mass of the clouds by using the depth of the events, duration, and 
standard extinction law. Thus, for an object like that in Fig. |3]F, with P~8 days and dip 
duration ~2 days, the size of the occulting body would be about 12 R^,. This large size is 
better associated with the warped disk model than with a discrete cloud in the disk. Given 
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the average J band depth of ~ 1 mag, for the narrowest dips (i.e., dip duration ~12 hr), 
the occuhing body must have a radius >0.6 R^,. If we consider the cl oud as a clump of d ust 



grains, we could estimate its mass using the N(H) to relation from lGagne et al.l (Il995[ ). If 
the cloud had a radius ~1 R^,, its mass would be ~ 2 x 10^^ g. If the cloud were as large as 5 
R*, then the mass of the cloud increases to ~ 2.5 x lO^'' g. These masses are small relative to 
the Moon (MassMoon=7.35xl0^^g), hence, these clouds are probably not the source material 
to form a planet. 

A problem with the disk obscuration interpretation arises if the flux dips are seen in 
diskless WTTs. Indeed, there are 3 ONC members not previously known to have IR ex- 
cesses (SOY 053524.15-044930.1, SOY 053529.02-050604.0, and SOY 053516.74-052019.7; 
see Table [H]) which have apparent AA-Tau type flux dips. The two first stars show just one 
apparent flux dip; the third shows 5 dips with no clear period. Examination of the SEDs 
of these stars shows that the latter shows evidence of having an IR excess at 5.8 and 8 /im 
(though not at 3.6 and 4.5 /im), consistent with it having a disk with a large inner hole. In 
this scenario, if the inner disk edge is sufficiently far from the star, there could be no excess 
at <5 /im, but there still could be material which could occult our view of the star and cause 
a flux dip. The other two stars could have even larger inner disk holes. 



5. Conclusions 

The YSOVAR Orion data provides the flrst large-scale survey of the photometric vari- 
ability of YSOs in the mid-IR. In Fall 2009, Spitzer/IRAC observed at 3.6 and 4.5 /xm a 
0.9 sq. deg. area centered on the Trapezium cluster twice a day for 40 consecutive days, 
producing light curves for 1249 known disked YSOs and ~ 820 young but generally disk- 
less WTTs, with typical accuracies ~3%. These data are publicly accessible through the 
YSOVAR website (http://ysovar.ipac. caltechiedu/l o'^ i^i Table |2l 

We flnd that 70% of the disked YSOs and 44% of the WTTs are variable. The WTT 
variations are mostly spot-like, while there is a wide range of variability types for the YSOs. 
Accordingly, 79% of the variable WTTs are periodic, while periods are detected for only 24% 
of the YSOs. On Table |4] we list 150 new periods found based on these data. We identify 
flve new PMS eclipsing binary candidates, including one already known spectroscopic binary, 
Theta Ori E, and another one that is fainter than all previously known ONC PMS eclipsing 
binaries. Coordinates and preliminary periods for the four new eclipsing binary candidates 
are given in Table [51 

We have discussed in detail a sample of the Orion variables that show short duration 
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flux dips (AA Tau events). The mean magnitudes of tliese stars plus tlie deptli of tlie dips 
and periods are given in Table |6l We interpret tliese events as the star being extincted by 
either clouds of relatively higher opacity in the disk atmosphere or geometric disk warps of 
relatively higher latitude which pass through the line of sight of the star. This hypothesis is 
consistent with both the duration and the wavelength dependent behavior of the events. 

We would like to thank the anonymous referee for her helpful comments. This work 
is based on observations made with the Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. 
This work is also based [in part] on observations obtained with WIRCam, a joint project of 
CFHT, Taiwan, Korea, Canada, France, and the Canada-France-Hawaii Telescope (CFHT) 
which is operated by the National Research Council (NRC) of Canada, the Institute National 
des Sciences de I'Univers of the Centre National de la Recherche Scientifique of France, and 
the University of Hawaii. 

Facilities: Spitzer (IRAC) 

REFERENCES 

Alencar, S. H. P., Teixeira, P. S., Guimaraes, M. M., et al. 2010, ArXiv e-prints 
Anathpindika, S. & Whitworth, A. P. 2008, A&A, 487, 605 
Attridge, J. M. & Herbst, W. 1992, ApJ, 398, L61 

Baraffe, I., Chabrier, C, Allard, F., & Hauschildt, P. H. 1998, A&A, 337, 403 
Bertout, C. 2000, A&A, 363, 984 

Bouvier, J., Grankin, K. N., Alencar, S. H. P., et al. 2003, A&A, 409, 169 
Cargile, P. A., Stassun, K. C, & Mathieu, R. D. 2008, ApJ, 674, 329 
Carpenter, J. M., Hillenbrand, L. A., & Skrutskie, M. F. 2001, AJ, 121, 3160 
Choi, P. I. & Herbst, W. 1996, AJ, 111, 283 
Fazio, G. G., Hora, J. L., Allen, L. E., et al. 2004, ApJS, 154, 10 

Fedele, D., van den Ancker, M. E., Petr-Gotzens, M. G., & RafaneUi, P. 2007, A&A, 472, 
207 



Furesz, G., Hartmann, L. W., Megeath, S. T., Szentgyorgyi, A. H., & Hamden, E. T. 2008, 
ApJ, 676, 1109 

Gagne, M., Caillault, J., & Stauffer, J. R. 1995, ApJ, 445, 280 
Getman, K. V., Feigelson, E. D., Grosso, N., et al. 2005, ApJS, 160, 353 
Gutermuth, R. A., Megeath, S. T., Myers, R C, et al. 2009, ApJS, 184, 18 
Haro, G. 1969, Boletin de los Observatorios Tonantzintla y Tacubaya, 5, 79 
Herbig, G. H. & Griffin, R. F. 2006, AJ, 132, 1763 
Herbig, G. H. & Kameswara Rao, N. 1972, ApJ, 174, 401 

Herbst, W., Bailer- Jones, C. A. L., Mundt, R., Meisenheimer, K., & Wackermann, R. 2002, 
A&A, 396, 513 

Herbst, W., Herbst, D. K., Grossman, E. J., & Weinstein, D. 1994, AJ, 108, 1906 
Herbst, W., LeDuc, K., Hamilton, C. M., et al. 2010, ArXiv e-prints 
Hillenbrand, L. A., Carpenter, J. M., & Skrutskie, M. F. 2001, ApJ, 547, L53 
Hillenbrand, L. A., Strom, S. E., Calvet, N., et al. 1998, AJ, 116, 1816 
Irwin, J., Aigrain, S., Hodgkin, S., et al. 2007, MNRAS, 380, 541 

Johnson, J. J., Gehrz, R. D., Jones, T. J., Hackwell, J. A., & Grasdalen, G. L. 1990, AJ, 
100, 518 

Jones, B. F. & Walker, M. F. 1988, AJ, 95, 1755 
Kovacs, G., Zucker, S., & Mazeh, T. 2002, A&A, 391, 369 
McNamara, B. J. 1976, AJ, 81, 845 

Parenago, P. P. 1954, Trudy Gosudarstvennogo Astronomicheskogo Instituta, 25, 1 

Parsamian, E. S., Chavira, E., & Gonzalez, G. 1993, Revista Mexicana de Astronomia y 
Astrofisica, 25, 71 

Plavchan, P., Gee, A. H., Stapelfeldt, K., & Becker, A. 2008a, ApJ, 684, L37 

Plavchan, P., Jura, M., Kirkpatrick, J. D., Cutri, R. M., & Gallagher, S. C. 2008b, ApJS, 
175, 191 



- 14 - 

Rayner, J., McLean, I., Aspin, C, & McCaughrean, M. 1989, MNRAS, 241, 469 
RebuU, L. M. 2001, AJ, 121, 1676 

Reipurth, B., Rodriguez, L. F., & Chini, R. 1999, AJ, 118, 983 

Scargle, J. D. 1982, ApJ, 263, 835 

Stassun, K. G., Mathieu, R. D., Mazeh, T., & Vrba, F. J. 1999, AJ, 117, 2941 

Stassun, K. G., Mathieu, R. D., & Valenti, J. A. 2006, Nature, 440, 311 

Stassun, K. G., Mathieu, R. D., Vaz, L. R R., Stroud, N., & Vrba, F. J. 2004, ApJS, 151, 
357 

Takahashi, S., Saito, M., Ohashi, N., et al. 2008, ApJ, 688, 344 

Tian, K. R, van Leeuwen, F., Zhao, J. L., & Su, C. G. 1996, A&AS, 118, 503 

Tobin, J. J., Hartmann, L., Furesz, G., Mateo, M., & Megeath, S. T. 2009, ApJ, 697, 1103 

Walker, M. F. 1978, ApJ, 224, 546 

Welch, D. L. & Stetson, P. B. 1993, AJ, 105, 1813 

Werner, M. W., Roellig, T. L., Low, F. J., et al. 2004, ApJS, 154, 1 

Yu, K. C., Bally, J., & Devine, D. 1997, ApJ, 485, L45+ 



This preprint was prepared with the A AS macros v5.2. 



- 15 - 



10.0 
10.2 
10.4 
^ 10.6 
^ 10.8 
11.0 
11.2 
11.4 

5.5 
5.6 
„ 5.7 
5.8 
5.9 
6.0 



a) Class I 
J053435.44-053959.1, 




t [3.6] IRAC 0[4.5] IRAC 



11.4 
11.6 
11.8 

12.0'° "T 
12.2 
12.4 
12.6 



10 20 30 
MJD-55127 



40 



? ° 

4 ° 



Class I, OMC-2 IRS1 
J053524.73-051 030.1 



• o * 



J [3.6) IRAC c|)[4.5] IRAC 



7.2 

cj 



7.4" 

7.5 

7.6 



10.2 
10.4 
10.6 
10.8 

11.0 
11.2 

8.0 
8.2 
, 8.4 
' 8.6 
8.8 
9.0 



b) Class II 

J05351 7.06-052858.2 



i [3.6) IRAC (t:[4.5] IRAC 



10.6 
10.8 



11.2 
11.4 



10 20 30 40 
MJD-55127 

^ o f ° o o f o. ^ 

« ° o o o ® 

' O 9(5 8 « ■ SO 

r ^" ^° °^ 'o 



10 20 30 
MJD-55127 



40 



Class II 
J053522.32-052414.2 

j .[3.6) |RAC (I)[4.5] IRAC , 

10 20 30 40 
MJD-55127 



8.6 



CD en 
9.0 CO ^ 



9.2 
9.4 



8.4 
8.6 
] 8.8 
9.0 
9.2 

11.2 
11.4 
11.6 
11.8 

12.0 
12.2 



?05351 6.33-052932.7 



j ,[3.6) |RAC (I>)4.5] IRAC 



10 20 30 
MJD-55127 



40 



f) 


Class II 

J053421 .39-04551 5.0 . 














i (3.6] IRAC 


J [4.5] IRAC 



9.4 

9.6 CD 

CO. 

9.8 
10.0 

12.0 
12.2 

CD 

12.42. 

12.6 

12.8 



10 20 30 
MJD-55127 



40 



Fig. 1. — Sample IRAC light curves for several stars showing the range of different variability 
types found in our data: a) slow drifts in brightness that may be caused by a slow change in 
accretion rate or by self-shadowing in a warped disk; b) rapid chaotic flux variations probably 
caused by changes in the accretion geometry; c) flares; d) periodic variations mixed with 
other longer term changes; e) periodic variations arising from photospheric spots and/or 
disk warps; f ) changes in color that may be produced by variable extinction. The symbols 
are •: [3.6]; o: [4.5]. [3.6] and [4.5] magnitudes have been plotted in the right and left vertical 
axis respectively. 
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Fig. 2. — Sample light curves for several stars showing shorter wavelength ground-based 
as well as Spitzer data selected to display the different variability types found in our data: 
a) larger amplitudes at shorter wavelengths; b) similar amplitudes of variation at all wave- 
lengths; c) larger IRAC amplitudes than the shorter wavelength amplitudes; d) IRAC varia- 
tions that are no n- apparent at shorter wavelengths; e) difference in phase between the optical 
an d near-IR data relat ive to IRAC data. Panel f) shows the phased light curve of star 1165 
of IStassun et al.l (119991 ). The symbols are •: [3.6]; o: [4.5]; *: J (dark green: UKIRT, light 
green: CFHT); and -|-: Ic (pink: SLOTIS, purple: APO). [3.6] and [4.5] magnitudes have 
been plotted in the right and left vertical axis respectively. Ic and J light curves have been 
shifted in the y-axis to match the mean IRAC values. The amount shifted is stated in each 
panel. 
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Fig. 3. — Sample light curves for six of the stars showing flux dips that can be discerned 
against a steady or slowly varying continuum. The two leftmost objects show broad dips 
that are not repeated over the duration of the observation. The two central panels show two 
of the narrowest dips and the two rightmost sources show periodic dips. The symbols are 
• : [3.6]; o: [4.5]; o: K,; *: J (dark green: UKIRT, light green: CFHT); and +: 1^ (pink: 
SLOTIS, purple: APO). [3.6] and [4.5] magnitudes have been plotted in the right and left 
vertical axis respectively. Ic, J, and Ks light curves have been shifted in the y-axis to match 
the mean IRAC values. The amount shifted is stated in each panel. 



Table 1. Summary of observations 



Telescope/Instrument 


Filter 


istart/iind dates 


N epochs 


Central coordinates (J 2000) 


Area (arcmin) 


Exptime/epoch (sec.) 


Spitzcr/IRAC 




Oct. 23 - Dec. 1, 2009 


81 


05;35;10 -05;42;25 


30x35 


(10.4&0.4)x4 


Spitzcr/IRAC 




Oct. 23 - Dec. 1, 2009 


81 


05;36;06 -05;16;22 


5x25 


(10.4&0.4)x4 


Spitzcr/IRAC 


[3.6], [4.5] 


Oct. 23 - Dec. 1, 2009 


81 


05;35;21 -04;49;47 


30x35 


(10.4&0.4)x4 


Si)itzcr/1RAC 




Oct. 23 - Dec. 1. 2009 


81 


05:.3'1:37 -05:20:00 


10 X 25 


(10.'ifr0.4)x4 


Spilzcr/IR AC 




Oct.. 2:i - Dec. 1, 2009 


81 


0.''):3.''):2i -()■''): 17:32 


20 X 2'i 


1.2x20 


UKIRT/WFCAM 


J 


Oct. 20 - Dec. 22, 2009 


32 


05:35:17 -05:22:49 


52x104 


2x3 


CFHT/WIRCam 




Oct. 27 - Nov. 8, 2009 


11 


05:35:14 -04:41:34 






CFHT/WIRCam 




Oct. 27 - Nov. 8, 2009 


11 


05:35:14 -05:02:14 






CFHT/WIRCam 


J,Ks 


Oct. 27 - Nov. 8, 2009 


11 


05:36:04 -05:22:46 


21x21 


(5x2)x7(J), 5x7(Ks) 


CFHT/WIRCam 




Oct. 27 - Nov. 8, 2009 


11 


05:34:26 -05:22:46 






CFHT/WIRCam 




Oct. 27 - Nov. 8, 2009 


11 


05:35:14 -05:43:28 






NMSU APOlm/CCD 


Ic 


Oct. 24, 2009 - Feb. 09, 2010 


27 


05:35:18 -05:24:03 


15x15 


60x20 


SuperLotis/CCD 


Ic 


Oct. 28 - Dec. 12, 2009 


23 


05:35:11 -05:41:02 


16x16 


180x12 



- 19 - 



Table 2. Time series at the [3.6] and [4.5] IRAC bands for the ONC stars. (Full table 
available in the electronic version of this article.) 





()l>j(-Ct 






AMD (days) 


IHAC baud 


].uag 




R\ .12000 (dc-R) 


Dec. .12000 (dc-g) 


ISOY. 


.J053532, 


.71- 


-045011, 


.9 


55127.50781250 


IRACl 


8.725 


0.006 


83.886268 


-4.836648 


ISOY. 


.J053532, 


.71- 


-045011, 


.9 


55127.76562500 


IRACl 


8.715 


0.006 


83.886268 


-4.836648 


ISOY. 


.J053532 


.71- 


-045011 


.9 


55128.00781250 


IRACl 


8.704 


0.006 


83.886268 


-4.836648 


ISOY. 


.J053532, 


.71- 


■045011, 


.9 


55128.37890625 


IRACl 


8.741 


0.004 


83.886268 


-4.836648 


ISOY. 


.J053532, 


.71- 


■045011 


.9 


55128.76171875 


IRACl 


8.814 


0.006 


83.886268 


-4.836648 


ISOY. 


.J053532, 


.71- 


■045011 


.9 


55129.25781250 


IRACl 


8.795 


0.005 


83.886268 


-4.836648 


ISOY. 


.J053532, 


.71- 


■045011, 


.9 


55129.89843750 


IRACl 


8.823 


0.006 


83.886268 


-4.836648 


ISOY. 


.J053532, 


.71- 


■045011, 


.9 


55130.69921875 


IRACl 


8.905 


0.005 


83.886268 


-4.836648 


ISOY. 


.J053532, 


.71- 


■045011, 


.9 


55131.35937500 


IRACl 


8.925 


0.007 


83.886268 


-4.836648 


ISOY. 


.J053532 


.71- 


■045011, 


.9 


55131.58593750 


IRACl 


8.912 


0.006 


83.886268 


-4.836648 
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Table 3. Statistics 





Class I 


Class 11 


Class 111 


# sources 


126 


1123 


820 


IRAC Variables 


106 


787 


366 


Variables w/ [3.6] amp. >0.2 mag 


71 


465 


92 


Variables w/ [3.6] amp. >0.5 mag 


15 


59 


8 


Variables w/ [3.6] amp. >1 mag 


3 


1 





[3.6] — [4.5] Color Variables w/ amp. >0.2 mag 


3 


5 





Prev. known Periods 


8 


266 


334 


Recovered Periods 


5 


125 


189 


Xew Periods 


8 


7() 


98 


# AA Tau analogs 


1 


37 


3 



Table 4. New periods for Orion members 



Source^ Period (days) p-value 



ISOY. 


.J053427.70-053155.4 


3.37 


4.78E-4 


ISOY. 


.,1053434.88-044243.5 


9.36 


1.31E-4 


ISOY. 


.J053435. 15-053210.4 


5.18 


2.46E-5 


ISOY. 


.J053435.98-045218.0 


2.34 


2.80E-3 


ISOY. 


.J053439. 11-053839.2 


2.41 


3.01E-2 


ISOY. 


.J053439.76-052425.6 


16.97 


1.65E-4 


ISOY. 


.J053441. 18-054611. 7 


8.02 


1.18E-3 


ISOY. 


.J053441. 71-052653.0 


8.47 


9.13E-2 


ISOY. 


.J053441. 87-055502.6 


1.66 


3.23E-6 


ISOY. 


.J053442. 74-052837.6 


6.15 


1.40E-7 


ISOY. 


.J053445.00-045559.2 


3.40 


3.60E-3 


ISOY. 


.J053447.64-054350.7 


4.27 


1.35E-2 


ISOY. 


.J053447.74-052632.1 


1.84 


5.73E-6 


ISOY. 


.J053447.85-044228.9 


8.36 


8.55E-4 


ISOY. 


..1053449.57-052903.4 


14.73 


9.90E-7 


ISOY. 


.J053450.66-050407.7 


9.24 


2.03E-3 


ISOY. 


.J053450.72-045836.8 


6.32 


3.80E-7 


ISOY. 


.J053452.22-045102.7 


5.15 


1.52E-2 


ISOY. 


.J053452.49-044940.4 


6.25 


5.78E-2 


ISOY. 


.J053452.60-051536.6 


2.35 


3.00E-7 


ISOY. 


.J053453.12-055955.0 


2.55 


3.09E-2 


ISOY. 


.J053453.99-044537.3 


2.38 


l.OOE-8 


ISOY. 


.J053454.83-052512.6 


3.75 


2.88E-4 


ISOY. 


.J053454.90-054643.9 


13.55 


4.70E-7 


ISOY. 


.J053454.92-050128.4 


13.90 


2.89E-2 


ISOY. 


.J053457.20-050823.9 


1.63 


2.85E-5 


ISOY. 


.J053458.53-060000.4 


5.98 


6.03E-2 


ISOY. 


.J053500.22-052409.2 


7.29 


2.60E-4 


ISOY. 


.J053500.54-054859.1 


6.63 


9.01E-6 


ISOY. 


.J053501.08-052304.1 


3.98 


6.,39E-5 


ISOY. 


.J053501. 10-052337.3 


7.12 


4.55E-3 


ISOY. 


.J053501.34-052811.9 


7.23 


7.00E-8 


ISOY. 


.,1053501.34-054113.5 


8.41 


1.36E-5 


ISOY. 


.,1053502.38-051548.0 


6.09 


1.49E-6 


ISOY. 


.J053503.12-050917.0 


4.32 


5.17E-5 


ISOY. 


.J053503.24-051726.4 


4.73 


5.13E-6 


ISOY. 


.J053503.70-052245.7 


5.11 


l.OOE-8 


ISOY. 


.J053504.56-052013.9 


8.43 


1.40E-7 


ISOY. 


.J053504.61-045829.0 


4.67 


1.04E-3 


ISOY. 


.J053504.80-060020.8 


5.13 


1.83E-3 


ISOY. 


.J053504.95-052109.2 


10.12 


2.00E-8 


ISOY. 


.J053505.45-052230.5 


3.56 


2.25E-4 


ISOY. 


.J053505.71-052354.1 


20.48 


5.41E-1 


ISOY. 


.J053505.79-053827.8 


2.80 


3.50E-7 


ISOY. 


.J053505.85-044843.3 


6.16 


1.63E-2 


ISOY. 


.J053506.77-055101.3 


6.76 


1.50E-4 
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Table 4 — Continued 



Source* Period (days) p-value 



ISOY. 


.J053506.86-051133.2 


5.47 


1.04E-6 


ISOY. 


.J053507.01-050134.8 


2.83 


8.27E-3 


ISOY. 


.J053507.70-055749.5 


3.75 


3.89E-2 


ISOY. 


.J053508.02-053244.3 


6.87 


2.13E-4 


ISOY. 


.J053508.07-054853.8 


7.97 


2.78E-4 


ISOY. 


.J053508.26-055000.3 


6.83 


1.04E-6 


ISOY. 


.J053509.95-051449.9 


2.19 


8.94E-6 


ISOY. 


.J053510.02-051816.6 


9.37 


l.OOE-8 


ISOY. 


.J053510.99-051521.8 


11.29 


3.53E-3 


ISOY. 


.J053511. 10-051601. 8 


4.75 


2.00E-8 


ISOY. 


.J053511.48-052352.0 


22.37 


1.49E-2 


ISOY. 


.J053511. 56-052448.0 


3.74 


7.68E-2 


ISOY. 


.J053511.62-051912.3 


6.77 


2.36E-2 


ISOY.J053511.89-052002.0 


4.66 


1.40E-6 


ISOY. 


.J053511.95-052845.0 


3.11 


1.13E-4 


ISOY. 


.J053513.46-053502.8 


7.28 


9.36E-4 


ISOY. 


.J053513.57-053508.1 


5.73 


3.24E-2 


ISOY. 


.J053514.66-050312.6 


7.66 


2.28E-4 


ISOY. 


.J053514.88-055142.2 


4.27 


6.33E-4 


ISOY. 


.J053515.48-052722.7 


6.54 


2.84E-6 


ISOY. 


.J053515.49-053511.9 


5.76 


3.71E-2 


ISOY. 


.J053515. 52-055316.2 


10.43 


6.00E-8 


ISOY. 


.J053515. 56-045308.1 


13.22 


3.90E-7 


ISOY. 


.J053515.65-055214.6 


8.67 


3.49E-4 


ISOY. 


.J053515.76-052309.9 


9.90 


3.00E-8 


ISOY. 


.J053516.33-052932.7 


16.79 


3.08E-6 


ISOY. 


.J053516.98-053246.4 


3.77 


1.42E-2 


ISOY. 


.J053517.10-051900.8 


17.46 


2.17E-3 


ISOY. 


.J053517.21-044113.6 


1.05 


6.44E-2 


ISOY. 


.J053517.35-052235.8 


7.10 


2.35E-2 


ISOY. 


.J053517.40-045957.2 


3.00 


9.60E-2 


ISOY. 


.J053517.53-051929.0 


6.72 


2.48E-5 


ISOY. 


.J053517.54-051613.1 


4.06 


7.10E-7 


ISOY. 


.J053517.97-054937.5 


6.27 


5.94E-2 


ISOY. 


.J053518.03-052205.4 


5.63 


1.06E-5 


ISOY. 


.J053518.23-051745.0 


3.28 


4.39E-6 


ISOY. 


.J053518.29-050805.0 


4.85 


2.71E-4 


ISOY. 


.J053518.70-051801.8 


7.37 


1.97E-4 


ISOY. 


.J053519.86-053103.8 


4.23 


6.85E-6 


ISOY. 


.J053520.02-052911.9 


2.79 


1.37E-2 


ISOY. 


.J053520.06-050358.5 


12.41 


9.96E-4 


ISOY. 


.J053520. 19-052308.6 


9.02 


3.32E-3 


ISOY. 


.J053520.29-054639.9 


7.90 


5.73E-2 


ISOY. 


.J053520. 72-051926.5 


22.21 


l.OOE-8 


ISOY. 


.J053521.38-050942.3 


3.76 


9.00E-8 


ISOY. 


.J053521.62-052325.8 


13.72 


3.39E-5 
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Table 4 — Continued 



Source* Period (days) p-value 



ISOY. 


.J053522. 10-051857.6 


5.78 


8.76E-6 


ISOY. 


.J053522. 32-041133. 3 


10.29 


9.60E-2 


ISOY. 


.J053522. 73-051838. 2 


5.75 


5.20E-7 


ISOY. 


.J053522. 98-051521.8 


4.32 


1.33E-3 


ISOY. 


.J053523.04-052941.5 


2.94 


1.30E-7 


ISOY. 


.J053523. 18-052228.3 


8.33 


5.18E-3 


ISOY. 


.J053523.97-055942.0 


18.70 


7.13E-5 


ISOY. 


.J053524.33-050120.5 


14.65 


1.23E-4 


ISOY. 


.J053524.34-052232.3 


18.91 


3.02E-6 


ISOY. 


.J053524.62-052104.2 


16.78 


7.10E-3 


ISOY. 


.J053524.67-044943.0 


17.15 


9.00E-7 


ISOY. 


.J053524.73-051030.1 


4.37 


3.47E-4 


ISOY. 


.J053525. 18-050509.3 


5.64 


9.42E-5 


ISOY.J053526.06-052121.0 


2.58 


1.28E-4 


ISOY. 


.J053526.47-053016.4 


3.70 


4.34E-5 


ISOY. 


.J053526.88-044730.7 


3.91 


3.07E-5 


ISOY. 


.J053527.00-054845.8 


11.27 


4.94E-2 


ISOY. 


.J053527.27-050527.0 


17.86 


6.35E-3 


ISOY. 


.J053527.66-054255.1 


12.31 


6.22E-4 


ISOY. 


.J053528.19-050341.2 


23.52 


3.48E-5 


ISOY. 


.J053529.65-052002.2 


2.44 


4.66E-6 


ISOY. 


.J053530.61-045936.1 


8.57 


3.05E-4 


ISOY. 


.J053532.67-054528.4 


2.53 


6.30E-7 


ISOY. 


.J053532.71-045011.9 


14.38 


4.32E-6 


ISOY. 


.J053532.96-051204.7 


1.30 


1.35E-4 


ISOY. 


.J053533.34-051145.6 


16.19 


1.40E-7 


ISOY. 


.J053533.60-052209.8 


21.48 


1.20E-7 


ISOY. 


.J053534.03-055505.8 


11.47 


1.04E-4 


ISOY. 


.J053535. 16-051946.9 


4.12 


1.12E-2 


ISOY. 


.J053535. 22-044739.6 


10.24 


3.88E-2 


ISOY. 


.J053536.49-044259.2 


2.46 


2.34E-2 


ISOY. 


.J053536.50-052009.4 


9.87 


3.97E-4 


ISOY. 


.J053538.56-050803.4 


16.77 


3.49E-6 


ISOY. 


.J053539.07-052025.7 


4.51 


1.82E-2 


ISOY. 


.J053539.49-044019.3 


3.31 


3.53E-3 


ISOY. 


.J053539.74-045141.7 


3.05 


6.27E-7 


ISOY. 


.J053539.77-044024.3 


17.91 


6.18E-3 


ISOY. 


.J053540.15-053723.7 


6.31 


4.05E-2 


ISOY. 


.J053541. 70-052014.9 


6.67 


1.72E-2 


ISOY. 


.J053541. 79-045027. 7 


24.80 


1.26E-3 


ISOY. 


.J053542.01-051011.5 


8.21 


6.91E-4 


ISOY. 


.J053544.09-050837.6 


4.46 


1.40E-7 


ISOY. 


.J053545.32-044334.9 


2.36 


4.33E-2 


ISOY. 


.J053546.23-051808.6 


9.71 


2.00E-8 


ISOY. 


.J053547.80-051030.8 


2.26 


2.47E-6 


ISOY. 


.J053551.07-051508.9 


9.36 


4.56E-5 



Table 4 — Continued 





Source^ 


Period (days) 


p-value ^ 


ISOY. 


.J053552. 10-044915.1 


3.76 


9.46E-4 


ISOY. 


.J053553.21-044734.7 


1.73 


2.99E-2 


ISOY. 


.J053554.63-052707.5 


8.45 


6.27E-7 


ISOY. 


.J053556.10-052533.0 


5.17 


2.46E-5 


ISOY. 


.J053558.94-053253.9 


8.31 


8.38E-4 


ISOY. 


.J053559.17-053552.7 


19.04 


1.54E-3 


ISOY. 


.J053605.95-050041.2 


3.56 


6,04E-3 


ISOY. 


.J053606.67-045512.1 


1.66 


l.OlE-4 


ISOY. 


.J053606.99-051334.2 


2.65 


5.66E-6 


ISOY. 


.J053609.28-045000.8 


15.16 


1.07E-2 


ISOY. 


.J053612.32-045819.0 


4.09 


9.61E-2 


ISOY. 


.J053614.94-051533.7 


4.72 


4.05E-3 


'"ISOY stands for Initial Spit 


zer Orion YSO. 





The p-value is the probability of chance occurrence re- 
turned by the periodogram service. 



Table 5. New PMS Eclipsing Binary candidates 



Source* Period (days) 



ISOY.J053526. 88-044730. 7 3.91 

ISOY_.J053518. 03-052205. 4 5.63 

ISOY_,T053505. 71-052354.1 20.39 

ISOY.J053605.95-050041.2 3.56 



''ISOY stands for Initial Spitzer Orion YSO. 



Table 6. ONC AA Tau analogs 



Source^ 


mean I 


mean J 


mean [3.6] 


mean [4.5] 


Dip depth 


Dip depth J*^ 


Dip depth [3.6 ]^ 


Dip depth [4.5]^ 


Dip Period 
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mag 
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ISOY_053510.94^043957.6 




13.636 


11.342 


11.075 




1.25 


0.3 


0.3 




ISOY_053510.99-051521.8 




16.255 


9.17 


8.561 




1.7 


0.6 


0.35 


ll.29ib0.56 


ISOY_053516.33-051538.0 




11.961 


9.617 


9.354 




0.6 


0.23 


0.20 




ISOY_053516.74-052020.0 




13.275 


9.621 


9.34 




0.8 


0.2 


0.18 




ISOY_053518.42-044000.2 




13.522 


12.018 


11.67 




0.28 


0.12 






ISOY_053518.70-051801.8 


12.457 


13.932 


10.153 


9.85 


0.7 


1.1 


0.4 


0.4 


7.37±0.80 


ISOY_053523.98-052509.8 


10.991 


13.227 


10.6 


9.774 


0.35 


0.35 


0.17 


0.11 




ISOY_053524. 14-044930.3 




13.941 


11.904 


11.515 




0.79 


0.46 


0.44 




ISOY_053525.51-054544.8 


14.306 


12.887 


9.816 


9.38 


1.6 


1.25 


0.30 


0.25 


7.57ib0.31 


ISOY_053526.73-051645.1 


11.22 


13.185 


11.006 


10.806 


1.1 


0.6 


0.26 


0.24 


2.96it0.62 


ISOY_053529.03-050604.1 




11.675 


10.797 


10.802 




0.24 


0.06 


0.06 





Table 6 — Continued 



Source* mean I mean J mean [3.6] mean [4.5] Dip depth Dip depth j'' Dip depth [3.6 J*' Dip depth [4.5]'= Dip Period 





mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


days 


ISOY.053530.45-052811.3 


12.139 


13.72 


11.575 


11.179 


0.46 


0.22 


0.12 


0.10 




ISOY_053534.41-051838.6 




15.817 


11.332 


10.781 




1.15 


0.50 


0.37 


4.75±0.33 


ISOY.053535. 61-053908. 1 


15.154 


13.793 


11.142 


10.653 


1.25 


1.15 


0.6 


0.5 




ISOY_053536. 42-050115. 5 




10.665 


8.632 


8.355 




0.69 


0.15 


0.14 




ISOY.053547.65-053738.8 




12.808 


10.413 


10.068 




0.95 


0.3 


0.25 




ISOY.053556. 10-052533.0 




15.613 


12.045 


11.727 




0.8 


0.35 


0.26 


5.17±0.15 


ISOY_053558.94-053253.9 




14.063 


11.466 


10.99 




1.05 


0.45 


0.42 





*ISOY stands for Initial Spitzer Orion YSO. 

''Dip depths are mean values for sources with more than one dip. 



